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Abstract 

We report on a lattice QCD calculation of the 1 = 2 tttt scattering length using the overlap 
fermion formulation for both sea and valence quarks. We investigate the consistency of the lattice 
data with the prediction of the next-to-next-to-leading order chiral perturbation theory after cor- 
recting hnite volume effects. The calculation is performed on gauge ensembles of two-flavor QCD 
generated by the JLQCD collaboration on a 16 3 x 32 lattice at a lattice spacing ~ 0.12 fm. 
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I. INTRODUCTION 



In Quantum Chromodynamics (QCD), the pattern of chiral symmetry breaking in its 
vacuum governs the interaction among pions, as dictated by the chiral effective theory jl]. 
Beyond the limit of massless soft pions, the pion interaction receives corrections due to finite 
pion mass and momentum, which can be systematically calculated using the framework of 
chiral perturbation theory (ChPT) |2|. For the tttt scattering, the two-loop calculation has 
been performed and the analysis of experimental data has been attempted so far {3]. 

Although the phenomenological analysis has provided evidence that the pion scatterings 
are fitted well by ChPT, it has its own limitation that the range of applicability of the chiral 
expansion is not known a priori. Since the quark mass in the nature is a constant, it is simply 
not possible to study the mass range that ChPT can be used. This is relevant to the question 
whether kaon can also be treated in ChPT on the equal footing. Furthermore, so-called the 
low-energy constants, the parameters appearing in ChPT, cannot be determined within 
ChPT. Eventually, one has to solve QCD in order to make a parameter-free comparison 
with the experimental data. 

Lattice QCD calculation provides a powerful tool to approach such a goal. As far as 
the convergence of the chiral expansion is concerned, the best studied quantities are the 
pion mass m n and decay constant f n (for a recent review, see for instance {4]). Among 
other lattice studies, the JLQCD and TWQCD collaborations investigated the convergence 
in two-flavor QCD, using the next-to-next-to-leading order (NNLO) of ChPT formula with 
various expansion parameters that differ only beyond next-to-leading orders (NLO) jlj]. 
The advantage of this particular work is in the use of the overlap fermion formulation on 
the lattice. Since the overlap fermion preserves exact chiral symmetry (in massless QCD) 
while maintaining the flavor symmetry as in the continuum theory, the application of the 
chiral effective theory is justified, in contrast to other lattice studies using the Wilson- type 
or staggered lattice fermions, for which some modification of ChPT with extra unknown 
constants is mandatory. The study |j] showed that the chiral expansion converges well up 
to the kaon mass region as long as the expansion parameter £ = m^/(47r/ T ) 2 is used where 
finite quark mass correction is included in m n and f n . With this expansion parameter, one 
could effectively resum a part of higher order terms. 

Extension of such study to other physical quantities is important in order to obtain a bet- 
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3, |8] . There have been a number 
- Il2j in the past, and 
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ter idea on the overall applicability of ChPT. The JLQCD and TWQCD collaborations have 
so far calculated and analyzed the electromagnetic and scalar form factors of pion 6| using 
the NNLO formula in ChPT. For these quantities to be consistent with the corresponding 
phenomenological analysis, it turned out that the NNLO terms have to be included. It is 
crucial to extend the test of ChPT beyond the simplest quantities studied so far. 

In this work, we study the I = 2 tttt scattering length. The lattice calculation of this 
quantity is easiest among other hadron scattering amplitudes, since the quark-flow diagrams 
are limited to those without pair-creation or annihilation between initial and final two-pion 
states. One uses the formula provided by Liischer to relate the two-pion state energy with 
the scattering length, or in general scattering phase shift 
of lattice calculations with and without quenched approximation 
recently more realistic calculations with light dynamical quarks have also been performed 
13l . |l4j . We carry out a similar calculation but using the overlap fermion formulation on 
the lattice, that has exact chiral symmetry at finite lattice spacings. 

Our calculation is carried out on a 16 3 x 32 lattice with lattice spacing a ~ 0.12 fm. 

n 

This two-flavor QCD ensemble is generated by the JLQCD and TWQCD collaborations [15l | 
using the overlap fermion formulation for sea quarks. The sea quark mass ranges from m s /6 
to m s with m s the physical strange quark mass. In order to maintain exact chiral symmetry, 
the global topological charge Q of the SU(3) gauge field is fixed to its initial value in the 
simulation, which is typically zero. This is irrelevant in the infinite volume limit as the 
global topological charge should not affect local physics, but induces finite volume effect 
that scales as 1/V at finite space-time volume V 21]. 

Since the lattice volume (1.9 fm) 3 is not as large as those used in previous studies, e.q. 
[131 . 1141 ] , finite volume effects have to be carefully investigated. The effect due to the fixed 
global topology can be corrected following the strategy outlined in 21] . The conventional 
finite volume effect due to pions wrapping around the lattice can also be estimated using 



ChPT, as calculated 



or pion mass and decay constant 22j, |23|]. In this work we use the 



formula developed in 24j for the tttt scattering. 



This paper is organized as follows. After describing the details of the simulation setup 
in the next section, we discuss the method to extract the two-pion state energy in our 
calculation in Section IIII1 In Section IIVI we show numerical results of our simulations and 
discuss the correction for finite volume effects. 
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II. LATTICE SETUP 



We use the overlap fermion formulation 16|, ll7] for both sea and valence quarks. Its 



Dirac operator for a massive quark is written as 

/ Tfl \ / Tfl \ 

D{m) = [m + —J + [m - —J j 5 sgn[H w (-m )], (1) 

where Hw{— mo) is the standard hermitian Wilson-Dirac operator with a large negative mass 
term. We choose —mo = —1.6 in the lattice unit. The physical quark mass is controlled by 
m. For the gauge sector, we use the Iwasaki gauge action |18|] at j3 = 2.30 together with extra 
(irrelevant) Wilson fermions that suppress the near-zero modes of Hy/{— m ) 19J, so that 
the fermionic determinant does not have singularity and the Hybrid Monte Carlo (HMC) 
simulation becomes feasible. 

The numerical simulation is performed on a 16 3 x 32 lattice with two flavors of dynamical 
fermions as one of the main projects of the JLQCD and TWQCD collaborations 15]. The 
attice spacing a is determined as 0.1184(3)(17)(12) fm from the Sommer scale r = 0.49 fm 
2o| . The sea and valence quark masses are set to 0.015, 0.025, 0.035, 0.050, 0.070, and 0.100 
in the lattice unit, which correspond to the mass range between m s /6 and m s with m s the 
physical strange quark mass. 

For each sea quark mass, we calculate the two-pion state energy on about 200 gauge 
configurations, each separated by 50 HMC trajectories. (The number of the configurations 
used in the analysis is 191, 193, 187, 193, 193 and 187 for the sea quark mass 0.015, 0.025, 
0.035, 0.050, 0.070 and 0.100, respectively.) The JLQCD collaboration has calculated and 
stored the lowest 50 pairs of eigenvalues and eigenvectors of the overlap-Dirac operator for 
their dynamical configurations. We utilize them to precondition the overlap solver, which 
makes the calculation faster by an order of magnitude [lsl ]. 

On these ensembles, the global topological charge Q is fixed to zero during the HMC sim- 
ulations in order to suppress the occurrence of the unphysical near-zero modes of Hy/{— m ). 
This makes the application of the overlap operator D(m) much faster. The effect of artifi- 
cially fixed topological charge on the physical quantities can be understood and indeed be 
corrected as a finite volume effect of 0(1/V) 21], which will be discussed later. 
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III. CORRELATION FUNCTIONS 



we need the 



In order to extract the txtx scattering length through Liischer's formula 
energies of tx and txtx system in a finite volume. The energy of a hadron state h is extracted 
from the correlation function Ch{t; t ), which describes the temporal propagation of the state 
from to to t. The source and sink operators are chosen such that they have some overlap 
with the state of interest. In our calculation, we consider C n (t;t ) and C T7r (t;to) defined by 

C n (t;t ) = -l^(0|7r(t,f)V(t )|0), (2) 

X 

cU*;*o) = -^^<°K 7r ^ f ) 7r ^^) t ^o)^o)|o>, (3) 

where Tx(t,x) is an interpolating field of the pion, a local pseudo-scalar density. The source 
operator W(to) corresponds to a wall source spread over the space. For the calculation of 
the scattering length, only the zero-momentum states are necessary; the correlators eq. (T2J) 
and eq. ((3]) are projected onto the zero-momentum states. 

The operators Tx(t,x) and W(t) are written in terms of quark fields u and d as 

n(t,x) = d(t, x)7 5 u(t, x), (4) 

W{t) = Y J d{t,x) 1 b u{t,y). (5) 

x,y 

The wall source is used on gauge configurations fixed to the Coulomb gauge. 

For the 1 = 2 scattering amplitude we can consider a scattering of two 7r + 's, so that no 
creation or annihilation of quarks occur between to and t. The Wick-contracted quark-line 
diagrams are shown in Figure [TJ In addition to the "direct" contribution where the two- 
pion states do not exchange valence quarks (left), there is a valence-quark exchange diagram, 
which is called "crossed" (right). 

Since the temporal extent of the lattice is finite, the correlation functions are not simply 
an exponential function, C(£;0) ~ exp(— m^t), decaying in time by a rate m n , but have 
a structure that reflects the boundary condition. For the periodic boundary condition in 
the temporal direction, the two-point function has the form exp(—m n t) + exp(— m n (T — t)), 
where the second term describes the state propagating backward in time. Here, T is the 
temporal extent of the lattice. 

For the case of the four-point function, each pion may propagate in forward or backward 
directions in time, so that there are four distinct contributions to the correlation function as 
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FIG. 1: The Wick-contracted quark-line diagrams for the four-point function of the 1 = 2 tttt 
scattering. Only the diagrams propagating forward in time are shown for simplicity. Each line 
represents a quark propagator < q{x)q{y) > where q is either u quark or d quark. One end of the 
line with a dot represents a source q{y) and the other end of the line with an arrow represents a 
sink q(x). The left and right diagrams are called "direct" and "crossed", respectively. 
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FIG. 2: Four "direct" Wick-contracted quark-line diagrams for the four-point function of the 1=2 
7T7T scattering, which exist due to the periodic boundary condition in the temporal direction. 



shown in Figure |2j Two of them contain the two-pion state, while the others contain only 
one pion at a given time-slice. Thus, the four-point function C nn (t; 0) may have a form 

C„(t;i ) ~ a (e- E ^ + e -^ T -*)) + p e -m«t e - m „(T-t) 

= 2acosh(-E nn (t-T/2)) + (3e- m * T , (6) 

where a and (3 are constants. E nn represents the two-pion state energy. The second term in 
eq. (jSJ) represents the effect from two pions propagating in opposite directions in time and is 
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called the wrap-around term since one pion wraps around the temporal direction. It amounts 
to a constant contribution. Although this effect is suppressed for large T as exp(— m^T), 
its contamination could become non-negligible in the middle of the lattice, since the two- 
pion state signal is also suppressed as exp(— E V7T T/2) ~ exp(— m n T). We therefore fit the 
two-pion correlator by eq. 

In order to identify the time-separations where the ground state dominates, we consider 
the following ratios: 

C^t + 1) 



cosh(-m^(t+l-T/2)) 
^ cosh (-m^t- T/2)) ' 1 } 

C n7T (t + 1) — C-K-xit) 

(8) 



C n7T (t) ~~ C^it — 1) 

cosh (-E nn (t + 1 - T/2)) - cosh (-E n7T (t - T/2)) 



cosh (-E n7r (t - T/2)) - cosh (-E n7T (t - 1 - T/2)) ' 

At large time separations, these ratios approach the form given in the second line of each 
equation, as the ground state dominates. From these ratios, we may extract the effec- 
tive mass m n (t) and effective two-pion energy E n7T {t) for each time slice. By inspecting 
the resulting effective energies, we are able to identify the region where the ground state 
dominates. 

In the actual simulation, we adopt a technique called low mode average (LMA), Q,S] 
for the correlation function. In this technique the quark propagator is separated into the 
eigenmodes of the overlap-Dirac operator. Then, a part of the correlation function, in which 
only low modes of the quark propagator participate, is averaged over the time of the source, 
as shown in Appendix El By taking an average we expect to gain statistics and have a more 
stable plateau in the correlation function of pions. 

Showing numerical results of the effective energy, here we discuss the impacts of the LMA 
and the wrap-around effects on the correlation functions. 

For the two-point function, it is found that the correlation function with LMA has smaller 
statistical errors than the one without LMA for all quark mass. This is consistent with 



previous results 



381 ] . For the four-point function, we show in Figure 3 the two-pion effective 



energy obtained from eq. (JE1) with and without LMA. The results show that LMA works well 
for small quark masses, m = 0.015 — 0.050, but it is not so effective for larger quark masses. 
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FIG. 3: Comparison of the two-pion effective energy E nn (t) with LMA (circles) and without LMA 
(crosses) at each quark mass. For all quark masses E nn (t) is extracted with the wrap-around term. 



Especially, for m = 0.100, the data show larger statistical errors with LMA. We therefore 
use the effective energy without LMA at the heaviest quark mass for the following analysis. 

In Figure 4, we compare the effective energy extracted from the four-point function with 
and without taking account of the wrap-around term. The latter is done by fitting the 
effective energy only with the hyperbolic cosine term, the first term of eq. (Q. On the one 
hand the effective energy with the wrap-around term shows plateaus, while the effective 
energy without the wrap-around term rapidly decreases as the time increases. On the 
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FIG. 4: Comparison of the effective two-pion energy E vn (t) with the wrap-around term (circles) 
and without the wrap-around term (crosses) at each quark mass. For ma = 0.100 E nn (t) is 
extracted without LMA while for other quark masses E^it) is extracted without LMA. 



other hand the effective energy with the wrap-around term has larger statistical errors near 
t ~ T/2. This can be understood as follows. As t approaches T/2, the hyperbolic cosine 
behaves like a constant. Due to this nearly degenerate behavior, the errors become large 
when one fits the correlation functions near t ~ T/2 to the sum of the hyperbolic cosine and 
the constant terms. 

As the quark mass becomes smaller, the effective energy without the wrap-around term 
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TABLE I: Summary of the quark mass, m q , the ground state energies, m n (L) and E n7T extracted 
from and two- and four-point functions and the statistics. We note m^L) in order to distinguish 
it from the pion mass including the finite volume effect, m n . 





' ' "7T \ J 


E 


statistics 


0.015 


0.1716(12) 


0.3697(38) 


191 


0.025 


0.21984(87) 


0.4565(22) 


193 


0.035 


0.26104(84) 


0.5345(19) 


187 


0.050 


0.31332(80) 


0.6407(25) 


193 


0.070 


0.37101(85) 


0.7547(32) 


193 


0.100 


0.4475(15) 


0.9033(31) 


187 



starts to decrease earlier in time. It is, therefore, absolutely necessary to use the effective 
energy with the wrap-around term for smaller quark masses, m = 0.015 — 0.050. For larger 
quark masses, m = 0.070 and 0.100, however, the situation is different. The time region near 
t ~ T/2 turns out unimportant due to large statistical errors when one extracts the ground 
state energy, E nn , from the effective energy with the wrap-around term. In the time region, 
which is important in determining E n7T , the effect of the wrap-around term is small. For 
m = 0.070 and 0.100, it turns out that E n7T determined from the effective energy without 
the wrap-around term in the range 6 < t < 9 coincide with the one with the wrap-around 
term. In later analysis we use the former. 

In summary, using LMA for the whole quark masses except m = 0.100 and including the 
wrap-around terms except m = 0.700 and 0.100, we fit the correlation functions in the time 
range 9 < t < 15 and obtain the ground state energies m 7T and E 7T7T . For m = 0.100, the 
effective energy is extracted without the LMA, and the time range is chosen to be 6 < t < 9 
for m = 0.070 and 0.100 so that the effective energy coincides with and without the wrap- 
around term. The results are shown in Table [I] together with the statistics, i.e. the number 
of configurations used in the calculation of two- and four-point functions. 
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IV. SCATTERING LENGTH 



Next, we explain how we extract the scattering length at the physical pion mass from the 
energies m n and E nn in a finite volume obtained in the previous section. 



A. Liischer's Formula 

Liischer's formula (3, 0] relates the energy of two hadrons in a box of size L with the 
scattering phase shift of two hadrons 5. For the two-pion case, the relation is 



1 nf kL \ 

kcoto = — -o — 
nL V 2tt / 



(9) 



k is the pion momentum in the center-of-mass system and is related to the nir energy, E nn , 
as E W7T = 2^Jk 2 + m*. E W7T is measured on the lattice in the center-of-mass system, obtained 
in the previous section. The function S is defined as follows 



q=2-Kn/L 



hm 

A^oo -n't 

\n\<A ' b 



k 2 L 2 



d 3 q 
- AtvA 



q 2 — k 2 



(10) 
(11) 



Near the nir threshold, the inverse functions of S can be expanded in terms of 1/L and it 
leads to the following expression 



(I 



1 I 71 " I 

1 + Ci— — + C 2 



n \ 2 , 



L 



0(L- 



(12) 



m-jrL 3 [ L 
where a nn is the irir scattering length defined by 

= lim k cot 5. (13) 

The numerical constants C\ and c 2 are —2.837297 and 6.375183, respectively Q. 



B. Finite Volume Effects 

Through Liischer's formula eq. (112j) . we can extract the scattering length a T7r from m n and 
-E 1 ,,-^ measured on the lattice. However, there are corrections neglected in deriving eq. ([12"]) . 
i.e. the finite volume effects. These corrections are small if the volume is large enough. But 
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the spatial size of our lattice L = 16 x 0.1184 fm ~ 1.9 fm, which is not large enough to 
neglect these corrections a priori. We therefore investigate the impact of finite volume effects 
in our data. We investigate the conventional effect due to pions wrapping around the lattice 
and also the effect due to the fixed global topology. 



1. Analytical Formula 



Since the periodic boundary condition is imposed in the spacial direction of the lattice, 
particles can wrap around the lattice and quantum fluctuations of the virtual particles, 
i.e. loop contributions, to the physical quantities are modified. The modifications induce 
corrections of 0(exp(— mL)) to physical quantities measured in the finite box of size L, 
where m is the mass of the virtual particle. The corrections are dominated by the lightest 
particle, the pion, and become more important as the pion mass becomes smaller in the 
simulation. These pion-loop corrections can be estimated using ChPT. 

We need to know the correction for the hit scattering length Aa nn as well as the correction 
for the pion mass A%. Am T is known at NNLO of ChPT Q while Aa wn is obtained as 
described below. 

From the ttw scattering length in a finite box a nn (L) the tttt scattering length in the 
infinite volume be obtained as 



T7 i' rJT- (X7. 



/r\ ,t\ i. A(fccot<5) ' 
1 + m^{L)a^{L) ■ hm 



(14) 



where A(/ccot<5) is the correction for A; cot 5 in a finite box, 



A(fccot5)(L) = A; cot 5 - k cot S(L). 



(15) 



The formula for A (A; cot 5) is calculated at the NLO of ChPT in ref. 24j, which is expanded 
in (m^L)" 1 in the low- momentum limit as 



lim 

k-t-0 



A(k cot 6){L) 



57T 







dml 



iAX(m n ) + 2iAJ exp (4ml) 



|n|/0 



\n\m^L 
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FIG. 5: A (k cot 5) /m^ vs. m n L. The solid line represents the exact formula eq. (116j) while the 
dashed and dotted lines represent the expanded formulae with first two and three terms in eq. (j!7p , 
respectively. 



AI and AJ cxp are denned with the modified Bessel and Struve functions as in ref. 



24), and 

c(n) is the multiplicity factor, i.e. how many times the discretized momentum with n = \n\ 
appear in the sum. In Figure |5] we show A{k cot S) (L)/m n as a function of m n L. In this 
figure we show exact A(k cot 5) (L)/m n , eq. (fT6|) . and also A (k cot 5) (L)/m n expanded in 
(m T L) _1 with first two and three terms, respectively, as in eq. (fTT|) . From this figure, it is 
seen that the expansion of A(k cot S) (L)/m n does not converge well in the region m n L < 3. 
Therefore we use eq. ffTBl) instead of the expanded form, eq. ffT7|) . when we estimate the finite 
volume effect of virtual pion loops. 

In addition to these pion-loop corrections, we also consider the finite volume effect due to 
the fixed global topology. This effect causes corrections of 0(1/V) to the Green functions 
in general and can be estimated once the topological susceptibility and the Q dependence 
of the physical quantity of interests are known 2l|, |3j|. The topological susceptibility has 
been calculated for the JLQCD gauge configurations recently. The Q dependence can be 
obtained through ChPT. At the leading order of ChPT, only the pion mass depends on 
Q. In the present analysis, we therefore include this effect on the pion mass as done in 



13 



Then, the effect on the scattering length is calculated at the NLO of ChPT as explained 
in Appendix [HI The correction for — 2m. K {L) due to the fixed global topology can be 
expressed as 

KK. - 2m A L)) = ^ [^^| {I (in ?f> + l) + W ] + 0(V^] . (18) 

Here / is the pion decay constant in the chiral limit, Xt is the topological susceptibility, 
l' a , which is defined in Appendix El is a combination of the low-energy constants l{ at the 
renormalization scale /i, V is the four- dimensional volume V = L 3 x T, and m n (0) is the pion 
mass at 9 = 0. We obtain m n (6 = 0) by combining the pion mass calculated by ourselves for 



5, 



15. 



M- For / 

si I'M 



the fixed topology and the correction estimated by the JLQCD collaboration 
and Xt, we use the known results calculated on the same gauge configurations 
is calculated from the phenomenological values of the low-energy constants j^]. 

2. Numerical Results 

To estimate the finite volume effect, we first calculate E 7T7T — 2m 7T (L) from the uncorrected 
m n (L) and E nn (Tabled]). Secondly, we include finite volume effect due to the fixed global 
topology, 5(E nn — 2m 7r (L)), eq. (TT8l) . From E n7T — Im^iL) with corrections due to the fixed 
global topology taken into account, the scattering length a„ is obtained through Luscher's 
formula, eq. f lT2|) . Finally, we take into account the finite volume effect due to the pion- loop 
corrections, A(k cot 5)(L) using eq. (fT6j) . and obtain m^a^. 

In the above procedure we use the corrected pion mass including finite volume effects on 
the r.h.s. of eq. (TTSI) . eq. (fT21) and eq. (fTBT) . The corrected pion mass, m n , is given as 

m 7r (L) 



(19) 



(1 + R m )(l + T m )' 

where R m and T m are the corrections from the pion loop and the fixed global topology, re 



spectively, obtained by the JLQCD collaboration [35|. We summarize each of the corrections 
in Table [XTJ In this table, all the quantities are given in the lattice unit, and the first and 
second errors in 1 + T m are due to the errors of Xt and the low-energy constants, respectively. 

In Table II III we summarize finite volume effects on E n7T — 2m n (L) and the scattering 
length. From this table, it is seen that the pion-loop corrections, A(k cot S)(L)/m 7T , increase 
as the quark mass becomes lighter, while the corrections due to the fixed global topology, 
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TABLE II: Summary of finite volume effects on the pion mass. R m and T m are the corrections from 
the pion loop and the fixed global topology, respectively, obtained by the JLQCD collaboration 
351 ] . is the corrected pion mass. 



Vfln 
H 


1+ Rm 


1 + Tm 

l Til 




0.015 


1.0226(34) 


0.9761(08)(16) 


0.1718(14) 


0.025 


1.0095(19) 


0.9858(14)(13) 


0.2209(11) 


0.035 


1.0048(12) 


0.9889(21)(13) 


0.2626(10) 


0.050 


1.00188(51) 


0.99357(91)(89) 


0.31480(91) 


0.070 


1.00075(23) 


0.99515(41)(81) 


0.37245(92) 


0.100 


1.000214(70) 


0.9943(07)(11) 


0.4501(16) 



5[E W7T — 2m w (L)), become important for heavier quark masses. However, it turns out from 
771^(1/) a 7r7r (L) and m^a^ that these finite volume effects are not so large, in fact as large as 
a few %, except for the largest quark mass where the effect is about 10 %. In Table III II we 
also show the chiral expansion parameter, £ = m^/(167r 2 /^), which will be used in the chiral 
extrapolation in the next subsection. In our notation, the pion decay constant, should 
correspond to the experimental value, 132 MeV, at the physical pion mass. This should 
be kept in mind when one compares our results with others' because in some references a 
different convention for the pion decay constant is adopted. 



C. Chiral Extrapolation 

In this subsection we extrapolate the scattering length to the physical pion mass from 
those at the lattice data points obtained in the previous subsection. The overlap formalism 



ism are 
the 



has exact chiral symmetry on the lattice and the results obtained with this forma 
expected to be consistent with ChPT even at finite lattice spacings. Following ref. 
scattering lengths for different values of the chiral expansion parameter, £, are fitted by the 
expression of the scattering length at the NNLO of ChPT, which is cast in the form 

= -tt£ jl + £ (J ln£ + /g)) + e (- y OnO 2 + S h£ + «g) } + 0{i% (20) 
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TABLE III: Summary of finite volume effects on E nn — 2171^ and the scattering length. 5(E nw — 
2m w (L)^ and A ( fc c°t<5)(L) arg ^ e corrections from the fixed global topology and the pion loop, 
respectively. a^^L) and a nn are the uncorrected and corrected tttt scattering length, respectively. 
The chiral expansion parameter £ is also shown. 



m q 


- 2rn % (L) 


5(E n7T - 2m n (L)) 


A(fccot5)(L) 


m Tr(L)a n7T (L) m^a^ 




0.015 


0.0265(38) 


0.00016(16) 


0.0780(14) 


-0.204(23) -0.200(23) 


0.0545(18) 


0.025 


0.0168(20) 


-0.00003(15) 


0.03056(51) 


-0.220(22) -0.221(22) 


0.0782(15) 


0.035 


0.0127(16) 


-0.00016(16) 


0.01415(25) 


-0.238(25) -0.242(25) 


0.1000(21) 


0.050 


0.0140(21) 


-0.00027(15) 


0.005535(97) 


-0.361(45) -0.370(46) 


0.1208(22) 


0.070 


0.0125(19) 


-0.00036(14) 


0.002015(38) 


-0.447(56) -0.461(57) 


0.1507(21) 


0.100 


0.00836(49) 


-0.00087(27) 


0.000532(15) 


-0.453(23) -0.497(28) 


0.1828(35) 



as is given in Appendix Here l^ls are combinations of the low-energy constants in ChPT 
at a quark-mass independent scale. /S-s are taken as fitting parameters below. We use the 
scattering lengths at the lightest four to six quark masses for the fit. At the NNLO of ChPT, 
are taken as the fitting parameters in eq. (1201). while Z$ is kept to be zero as 
done in This is because practically l^} is hard to be determined. Instead, we checked 

(2) 

how much the fitting parameters change by including as its phenomenological value and 
incorporated the changes as systematic errors. For comparison, the data are also fitted by 
the NLO formula of ChPT with one parameter l^l by truncating 0(£ 3 ) terms in eq. ( !20l) . 

Figure 6 shows the results of the fit of the m^a^. The fits at the NLO (left hand side) 
and the NNLO (right hand side) of ChPT for the six quark masses are shown. The leading 
order (LO) result of ChPT is also shown there. At a glance we cannot conclude which of 
the LO, NLO and NNLO of ChPT fits the calculated m n a 7T7T best. Table IV summarizes 
the low-energy constants, l^} and l^, and m n a nn at the physical pion and the values of x 2 
obtained by the fit of ChPT. The values of x 2 are similar for the LO, NLO and NNLO fits 
when the results of lowest four or five quark masses are used. When the results of six quark 
masses are used, however, the NLO and NNLO fits have similar x 2 values but the LO fit 
has larger x 2 than the NLO and NNLO fits. This implies that the effect of the NLO and 
NNLO of ChPT appears only at the largest quark mass, m = 0.100. 



1(3 



TABLE IV: Summary of the results of the chiral extrapolation. 



order of ChPT 


data pts 


f'TTTV 


,(3) 


X 2 /d.o.f 


m-KCi-K-K (Physical) 


LO 


6 






3.1 


-0.04539 




5 






2.2 


-0.04539 




4 






2.8 


-0.04539 


NLO 


6 


1.97(23) 




1.7 


-0.04251(15) 




5 


2.39(44) 




1.8 


-0.04279(29) 




4 


2.18(60) 




2.5 


-0.04256(40) 


NNLO 


6 


5.8(1.2) 


-6.8(6.9) 


1.7 


-0.04410(69) 




5 


5.0(1.9) 


1(16) 


2.2 


-0.0437(11) 




4 


7.7(4.0) 


-27(38) 


3.2 


-0.0451(23) 



In Figure 7 we show Vkl and m^a^ at the physical pion mass, which are obtained by the 
fit at the NNLO and the NLO of ChPT. In the figure phenomenological values and the results 
of lattice simulations by other groups are also shown for comparison. The NNLO fit has 
larger errors in l^} and m^a^ than the NLO fit. This is because the former has more fitting 
parameters than the latter. From Figure [7] one sees that the central values of l^J obtained 
by the NNLO fit are consistent with the phenomenological value given in jsj], while, those 
obtained by the NLO fit seem to slightly deviate from the phenomenological values. These 
results including a change in convexity of the fitting curve suggest a significance of the NNLO 

(3) 

contributions especially for the heavier quark masses. Still, the values of the parameter, IW, 
which are determined using 4, 5 and 6 quark masses, are very much different. More precise 

(3) 

determination of demands more data points at lighter quark masses and careful analysis 
for the chiral extrapolation of the scattering length. It is noted, however, that the m^a^ 
can be extrapolated to the physical pion mass with relatively small errors. The extrapolated 
scattering length is compared with the newly reported experimental value extracted from 



the kaon decay measurement 



33 



34J. It is remarkable that the extrapolated value at the 



NNLO of ChPT agrees with the result of ref. 



We should also mention that our results 



are consistent with the results of other groups within errors. 
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FIG. 6: ChPT fit at NLO (left figure) and NNLO (right figure) level of the scattering length for 
the six quark masses. The x-axis is the expansion parameter £, and the y-axis is the scattering 
amplitude m n a nn . The straight line from the origin to the lower right is the leading order (LO) 
result of ChPT. In the left (right) figure the curved line shows the NLO (NNLO) fit of ChPT and 
the shaded area shows the la band of the fit. The upper left point is m^a^ extrapolated to the 
physical pion mass with using 6 quark masses. 
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FIG. 7: In} (left figure ) and 111^^^ at the physical scale (right figure). Crosses are our results 
obtained by the ChPT fit. The right (left) three points are the results of the NLO (NNLO) fit. 
4pt, 5pt, and 6pt below the points mean that 4, 5, and 6 quark masses are used for the fitting, 
respectively. The open circles are phenomenological values obtained by E865 [36], NA48/2 [37] and 
CGL [3] collaborations and the open triangles are the results of lattice simulations by NPLQCD 
group [13] and DESY group [14]. 
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V. SUMMARY 



We have calculated the I = 2 htx scattering length using the gauge configurations gen- 
erated by the JLQCD collaboration with the two-flavor dynamical overlap fermion. The 
overlap fermion action has exact chiral symmetry on the finite lattice, which enabled us to 
compare the calculated results with the chiral perturbation theory in the continuum. 

We have adopted the technique called low mode averaging in order to gain statistics for the 
correlation function and have taken into account the wrap-around term in the parametriza- 
tion of the correlation function in order to identify the contribution of the ground state in 
the correlation function. 

We have also investigated finite volume effects on the physical quantities, which appear 
due to pion-loop corrections and also due to fixed global topology. We have taken into 
account corrections of these effects on the ttw scattering length. In order to estimate the 
correction due to fixed global topology we have used the topological susceptibility obtained 
by JLQCD for the same gauge configurations. From our calculation, it is found that the 
corrections are small and the lattice volume used for the present calculation is large enough 
to suppress finite volume effects. 

The calculated results are fitted by the NNLO ChPT and also the NLO ChPT for com- 
parison, from which the scattering length at the physical pion mass is extrapolated. It is 
found that the NNLO ChPT and the NLO ChPT fit the calculated pion mass, the irw scat- 
tering length and the low-energy constant of ChPT equally well. However, the obtained 
scattering length and the low-energy constant of ChPT at the NNLO ChPT agree slightly 
better with the experimental ones than those obtained at the NLO ChPT. 

The scattering length extrapolated to the physical pion mass assuming the NNLO ChPT 
is m n a n7r = — 0.04410(69)(18), in which the numbers in the first and second parentheses 
represent the statistical error and the systematic error from using the truncated formula. 
The extrapolated scattering length is in good agreement with the experimental result. A 
combination of the low-energy constants of ChPT obtained by fitting the calculated results 
is also consistent with the phenomenological value. 
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Appendix A: Low Mode Average 

We decompose the quark propagator into the low mode ( "L" ) and the high mode ( "H" ) 
of the overlap-Dirac operator for each gauge configuration. Then, for example, the meson 
correlator is decomposed into four terms according to the two modes of the quark propagator 

as 

Or(t; t ) = C w (t; t ) + C HL (t; t ) + C LU (t; t ) + C LL (t; to). (Al) 

Among these four terms, Cll^^o) is the most important for physical quantities at low 
energies. The idea of LMA is to average only Cll(£; to) over the time of the source t 

1 T ~ l 

C<(t; t ) = Or(t; to) - C LL (t- 1 ) + - ^ C LL (t - * + *i! h). (A2) 

ti=0 

Eigenvalues and eigenvectors of the low modes of the overlap Dirac operator are stored and 
are used in order to calculate the average. We use the same technique for the two-pion state. 



Appendix B: finite volume effects from fixed global topology 



The n-point function with the topological charge Q fixed, is expressed as 



g2 = GJe a ) + 



1 d 2 G n {9) 



2 Xt V d6 2 

where G n {6) is the n-point function in the 6 vacuum, 



+ 



(Bl) 



G n {9) = A{6)e 



-m n (6)t 



(B2) 
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where Xt is the topological susceptibility, V is the volume of the lattice and 8 S = is the 
saddle point. 



Assuming that 



m n (0) = m n (O) + i^^ „-0 2 + --", (B3) 



(9=0 



we can approximate as 

G% ~ A Q e - m ™*, (B4) 

where 

< = ^(0) + f 1 - f-) + 0(V- 2 ) + • ■ ■ . (B5) 

Therefore, on the one hand E V1T — 2m n is given as 
Eg - 2m® = mi - 2m® (B6) 

= E „ {0) - 2nt ,(o) + ^(i-f) ^ {E " m d - 2mM) I +0 ,v- 



-2\ 



(B7) 

with the fixed topological charge, Q. On the other hand, E vn — 2m n is given from Liischer's 
formula, eq. (fT2|) . and the ChPT, eq. (1C3|) . as 

^ - 2m. = (1 + ©(IT 1 )) (B8) 

fl + T7 ! 4^f-^ + 0) + ---- (B9) 



4/^3 [ levrV. 2 V 2 



4/ 2 L 3 [ 16tt 2 / 2 V2 A* 2 
where 

'" = -^-5' 2 + 5''-5- (B11) 

In the last line, the following formula for the pion decay constant /„- is used. 

f^f^-Sni- 1 - 1 ')}*-' (B12) 



21 



where / is the pion decay constant in the chiral limit. Then, we assume that the 9- 
dependence of E n „ — 2m n appears only through that of the pion mass, which is taken 
into account by changing m n into the following form of m n (9) 



m 



n ->■ m n (9) = 2B m q cos(9/N f ), 



(B13) 



where m q is the quark mass and Nf is the number of the flavor. By taking the second 
derivative of eq. (IBlOj) by 9, substituting it into eq. (1B7|) and setting Q = 0, we obtain the 
following expression for the finite volume correction due to the fixed global topology: 



5(E nw - 2m n ) = (Eg - 2m£) - (^(0) - 2m ff (0)) 



Af 2 L 3 



1 ml(0) 



/ ln m|(0) 



(B14) 
(B15) 



Appendix C: Scattering length in ChPT at NNLO 



The 7T7T scattering length is given at the NNLO of ChPT in ref. 23] • The S-wave scattering 
length in / = 2 channel is expressed as 

"262 22tt 2 



tt£ {I - £[2 + h + 166 4 ] + e 



h646 4 | | + 0(( 4 ), (CI) 
to parametrize 



9 9 

where fej's are coefficients (multiplied by 167T 2 ) introduced in ref. 
the pion scattering amplitude, fej's are 

6 4 8~ 1 6 31 /4~ ~ ~ 47 

6i + I664 = -L + -/1 + 7J2 - -k ~ 2k - - + -tl 2 + & -k + 8/ 2 - h + 2/4 - - 
16 6 1 2 \ 3 12 



Vl6~~ 32- ~ 
+ £ irhk + —hh 



-r - 3Z 3 / 4 - 5/ 4 2 - ^/ 3 2 + 4/i + y Z 2 - ^3 - 6/4 + + fl + 16f 4 

(C2) 



with the low-energy constants Zj, fj in ChPT of chiral order 4 and 6 renormalized at scale \x 
and 



ml 



167TV 2 ' m /i 2 



Substituting eq. (IC2I) into eq. fjClj) and rearranging the result in the order of £ lead to 
the following expression: 



--L + lA + £ 2 (-y^ 2 + l bL + / c 



(C3) 
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with 

4~ 8~ l~ ~ 1 

~~ ~~ 3 ~~ 3 + 2 + ~~ 2' 

4~ ~ ~ 119 
h = — g/i— 8/2+^3-2/4 + — , 

l r , /16~ 32- ~ r \ r 4~ 16- 7 r ~ 163 22 „ 

c = 2 3 v T + T ~~ ~~ 5 / + 3 + T + 1 ~ + Iff ~ Y 7 " 1 ~~ 16f4 ' 

(C4) 

The scale dependence of fj is fixed so that 6j has no scale dependence, \idbijd\i = 0. There- 
fore, the right hand side of eq. (ICip . is scale independent, though truncated at the third 
order of £. Then, the right hand side of eq. (103|) as a whole is also scale independent. Thus, 
one can choose \x arbitrarily. When we fit the scattering length obtained from lattice simu- 
lations as a function of £, we want to make the fitting parameters quark-mass independent. 
Therefore, we choose \i = 4irf with the pion decay constant in the chiral limit / At this 
scale, it should be noted that 

L(^ = 47r/) = -ln£-ln|| 

= -ln£-2£/ 4 + 0(a 

= - hie - 2£/ 4 (/i = 4tt/) + 2£ln£ + 0(£ 2 ), (C5) 

in which we have used the chiral expansion of the pion decay constant f n = f{l+^U + 0(^ 2 )} 
with I4 being the low-energy constant at the pion-mass scale [si, 0] • 

Substituting eq. (1C5I) into eq. (1C3I) . setting the scale /i = 47r/ and sorting the result in 
order of £, we obtain the following expression for the pion scattering length M 

= |l + £ Q lne + + £ 2 (~f ( ln 0* - (4 + 3) hi£ + (/ c + 3Z 4 )) } + 0(£ 4 )- 

(C6) 

Note that all the low-energy constants are at the scale of fi = 4nf and are independent of 
the quark mass, so that we can take them as fitting parameters in the chiral extrapolation 
of the scattering length. 
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